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文章 编号 :1005-3085(2010)05-0771-10 


响应 变量 随机 缺失 下 的 变 系数 部 分 
线性 模型 的 经 验 似 然 推断 * 


赵 培 信 !， 薛 留 根 ? 
(1- 河池 学 院 数学 系 ， 广 西 宜 州 546300; 2- 北京 工业 大 学 应 用 数理 学 院 ， 北 京 100124) 
摘 要 : 本 文 考虑 了 响应 变量 随机 缺失 下 的 变 系 数 部 分 线性 模型 的 估计 问题 。 利 用 经 验 似 然 方法 ， 给 出 了 
参数 部 分 的 调整 经 验 似 然 比 函数 ， 证 明 其 渐 近 服从 标准 卡 方 分 布 。 进 而 构造 了 参数 部 分 的 置信 
域 ， 得 到 了 其 极 大 经 验 似 然 估计 的 最 优 参数 收敛 速度 和 渐 近 半 参 数 有 效 界 。 模 拟 结果 表明 调整 经 
验 似 然 方法 优 于 未 调整 的 经 验 似 然 方法 。 
关键 词 : 变 系数 部 分 线性 模型 ， 经 验 似 然 ， 置 信 域 ， 缺失 数 据 
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1 引言 


设 有 nn 个 个 体 的 样本 ，(Yi, Xi, Zi Ui), i=1, ,n， 本 文 考虑 如 下 的 变 系数 部 分 线性 模型 
Y; = X70(U;) + ZTB + éi, (1) 


其 中 Bb = (beba) Aq HARRIE, 0O) = (91(-),… CT A p HEAR Pl I. Yi X 
响应 变量 ，X; 为 p 维 协 变量 ，2 为 g 维 协 变量 ，Ui 为 一 维 协 变 量 。 不 失 一 般 性 ， 设 U0 在 区 
间 [0,1 上 取 值 。 并 且 假 定 B(ei | Ui, Xi, Zi) = 0。 

模型 (1) 包含 了 几 种 被 广泛 研究 的 模型 ， 例 如 线性 模型 83， 部 分 线性 模型 8-6 以 及 变 系数 
模型 (78) 等 都 是 模型 (1) 的 特殊 情况 。 对 模型 (1)， 在 数据 不 含有 缺失 的 情况 下 已 有 大 量 的 文献 
进行 研究 820， 但 是 在 实际 应 用 中 ， 响 应 变量 往往 因 某 种 原因 而 产生 缺失 。 当 响应 变量 含有 缺 
失 时 ， 现 有 的 统计 推断 方法 将 不 能 直接 应 用 。 如 果 仅 仅 用 可 以 完全 观测 到 的 样本 进行 统计 推 
断 ， 得 到 的 估计 往往 会 产生 偏差 ， 而 且 一 般 不 是 渐 近 有 效 的 。 关 于 缺失 数据 问题 的 研究 ， 目 前 
已 有 大 量 的 文献 对 缺失 数据 下 的 线性 模型 3 以 及 部 分 线性 模型 2,13] 进行 了 研究 。 但 是 关于 缺 
失 数据 下 的 变 系数 部 分 线性 模型 的 估计 问题 ， 目 前 还 没有 相关 的 文献 进行 研究 。 

本 文 在 响应 变量 随机 缺失 下 ， 利 用 经 验 似 然 方法 对 模型 (1) 参数 部 分 的 估计 问题 进行 研究 。 
给 出 了 参数 部 分 的 调整 经 验 似 然 比 函数 ， 并 证 其 渐 近 服从 标准 卡 方 分 布 ， 进 而 构造 了 O 的 置信 
域 。 另 外 本 文 还 证 明了 参数 部 分 的 极 大 经 验 似 然 估计 达到 了 最 优 的 参数 收敛 速度 并 且 是 渐 近 有 
效 的 。 最 后 通过 数据 模拟 ， 对 调整 的 经 验 似 然 方法 与 未 调整 的 经 验 似 然 方 法 的 有 限 样本 性 质 进 
行 了 比较 。 模 拟 结果 显示 本 文 的 方法 优 于 未 调整 的 经 验 似 然 方法 。 
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2 参数 部 分 的 调整 经 验 似 然 推 断 


BW (Yi, Xi, 2i,6i,0i), i51, ,n， 为 来 自 模型 (1) 的 一 个 不 完全 随机 样本 ， 其 中 XX;, 2Z; 和 Ui 
可 以 完全 观测 ， 当 6; = 1 时 ，Y; 可 以 观测 ，6; = 0 时 ， 玫 缺失 。 本 文 假定 为 随机 缺失 ， 即 
P(6=1|Y,X,2,U0)= P(6=1|X,2,U). 
由 模型 (1) 可知 
iY; = 6X7 O(Ui) + 5:27 B + ei, (2) 


其 中 ei = die, 为 独立 同 分 布 的 随机 变量 ， 并 且 Ele | Xi, Zi Ui) = 0. MH (2) 为 基于 完全 观 
测 数据 的 变 系数 部 分 线性 模型 ， 因 此 利用 文献 [9] 的 方法 ， 对 给 定 的 6， 最 小 化 (3) 式 则 得 
F O(u) 的 局 部 最 小 二 乘 估计 0b(u)。 

n p 

5 { Y; — ZTB — >> [ax + bx (Ui — u)] Xix } GiKnlu — Ui), (3) 

i=1 k=1 
其 中 Kal) =A K(-/h), K(-) ABR, HAH, Xn AX 的 第 天 个 元 素 。 记 

Z= (Zn, Za)”, Y = (Yi, Ya), 

五 为 p 阶 单位 阵 ，0;, 为 p 阶 0 矩阵 


S(u) = (Ip, Op) (D7 Uu Du) D] Qu 


Il 


(Si(u), ,Sn(u)), 
其 中 

Qu = diag( ô Kn (u = U1), KAR ,On Kp(u S Un)) 
为 nxn 对 角 阵 ， 并 且 


T 
Xi ds Xn 
Dy = 
( hi(U Re u)Xı LES ht, = u)Xn ) 
An x 2p 的 矩阵 ， 则 有 


Ô(U:) = >》 Se(Ui) (Ys — Ze P). (4) 
k=1 


把 (4) KARA (1) 式 ， 并 简单 计算 得 
Y; — X79(U,) ~ [Z; -— AU)" Xi] B + 4, (5) 
其 中 
flu) = > Sk(u)28，6(u) = >》 Se(u)¥e- 
k=1 


k=1 


引入 辅助 随机 向 量 
ði 





市 (8) = (Zi — (Ui)T Xi) {Yi — XP9(Ui) — (Zi — Ai)? Xi)" B}, (6) 


a 
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其 中 调整 因子 A(w) 为 A(u) = P(6; =1|U; = wv) 的 核 估计 ， 即 
A(u) = 》 wni(u)6, 
i=l 


这 里 
ne _ K((u-Ui)/h) 
i Ez- K((u — U;)/h) 


由 (5) 式 可 得 EAO} = o(1)。 根 据 这 一 信息 ，B 的 经 验 似 然 比 函数 可 定义 为 
R(8) = -2max{ > log(np,) |p. 20, Sop =1, 5 phe) = o}, 
i=l i=1 i=1 


其 中 mi = p(B), i= 1,… ,n。 对 任意 给 定 的 8， 假 设 0 在 点 (加 (8),… ,加 (6)) 所 构成 的 凸 集 
内 部 ， 则 RCO) 存在 唯一 的 解 。 利 用 Lagrange RAYE, TAH ÊO 表示 为 


R(B) = 2 log (1+ A7H:(B)), (7) 


t=1 


其 中 入 = A) 为 9 维 向 量 ， 且 满足 


M(B) 
LTH ® 
定理 1 设 第 4 节 中 的 条 件 C1-C5 成 立 ， 如 果 6 为 参数 真 值 ， 那 么 
R) = x2, 


其 中 A, RAMA ATRL, 2 为 自由 度 为 g 的 卡 方 分 布 。 
以 X2(1 一 @) 记 又 的 1 一 Qa 分 位 数 ，0 < a < 1。 由 定理 1 可 以 得 到 8 的 近似 1 一 a 置信 域 为 


Ca(B) = { B| R(B) < x2(1 — a) }. 
通过 最 小 化 (7) 式 可 以 得 到 68 的 一 个 估计 8, RO WERAKB BWA. jd 
nS a . T 
P= = 2 (Zi — Wi)" Xi) (Zi — Wi)" Xi), 


如 果 和 矩阵 工 是 可 着 的 ， 那 么 利用 Qin A Lawless!4) 的 定理 1 的 类 似 证 明 可 知 ，6 的 极 大 经 验 似 
然 估计 可 以 表示 为 


Ê = -1E SO (Z: — RO)T Xi) (Y: — XP 4) + op(n-4). (9) 
i=l 


令 
U() = E(GXP|U;=-), 06) = E(X:Z7 |U: =-), 0() = B(2|U; =), 


下 面 的 定理 给 出 了 6 的 浙 近 正太 性 。 


774 I E 数 学 学 R 第 27 卷 








定理 2 设 第 4 节 中 的 条 件 C1-C5 成 立 ， 那 么 


其 中 
X = E{ o?(U)[E(ZZ7) — ®(U)T LU) EV) }. 


由 定理 2 知 = 8 + Op(n-)， 即 AB 达到 了 最 优 的 参数 收 人 速度。 另外， 由 文献 [1 引 可 
知 马 为 半 参 数 有 效 信息 界 。 因 此 ， 本 文 提 出 的 调整 经 验 似 然 方法 所 得 的 估计 是 半 参 数 有 效 的 。 


3 ”数字 模拟 分 析 


为 实施 模拟 ， 我 们 从 模型 Y = sin(0.5xrU)X+1.52 十 e 中 产生 数据 ， 其 中 X ~ N(0,1) UA 
2Z 均 为 [2,2] 上 的 均匀 分 布 。Y 由 模型 产生 ， 其 中 模型 误差 e 服 从 均值 为 0， 方 差 为 0.5 的 正 态 
分 布 。 在 模拟 过 程 中 ， 取 n = 100， 实 验 重复 1000 次 。 并 且 缺 失 概率 Au) 分别 取 如 下 三 种 情 
况 来 代表 响应 变量 的 不 同 缺 失 水 平 : 

1) WRU) < 0.5, W A(u)=0.9+0.2U, FH A(u) = 0.95; 

2) 如 果 |U| < 0.5, WA(u) =0.7+0.2U, AM A(u) = 0.75; 

3) ”如果 |U| < 0.5, WA) =0440.2U, AM A(u) = 0.45- 

这 三 种 情况 对 应 的 数据 Y 的 缺失 概率 分 别 平均 为 10%，30% 和 60%。 核 函数 取 K(wu) = 
0.75(1 — 222)+， 用 “去 一 个 体 ” 交 叉 证 实 法 选取 带宽 hoy, HAE 


OVR) = SDa Y: -XFU - ZTA Y 
i=l 


达到 最 小 ， 其 中 6_;(Ui) A pa DH WER i PE O(U,) 和 6 的 估计 量 ，6_;(Ui) 可 通过 
Ht B_ RA (4) 式 得 到 。 

对 于 参数 6， 我 们 对 本 文 提出 的 调整 经 验 似 然 (ADEL) 方法 与 未 调整 经 验 似 然 方 法 (NAEL) 
进行 比较 。 后 者 是 只 用 完全 观测 到 的 样本 ， 直 接 把 (6) 式 中 的 调整 因子 A(U)-36; 用 6; 所 代 
蔡 ， 来 构造 经 验 似 然 比 函数 。 我 们 分 别 计算 了 1000 次 模拟 的 绝对 偏差 |6 一 6| 的 平均 值 (Bias), 
标准 误差 的 平均 值 (SE)， 和 置信 水 平 为 95% 的 置信 区 间 长 度 的 平均 值 (Len) 和 对 应 的 覆盖 概 
率 (Cov) 的 平均 值 。 模 拟 结果 见 表 1。 


Rl: 响应 变量 和 在 三 种 缺失 情况 下 的 模拟 结果 


ADEL NAEL 
1— E(A) Bias SE Len Cov Bias SE Len Cov 
10% 0.055 0.041 0.221 0.954 0.054 0.037 0.223 0.958 
30% 0.071 0.037 0.261 0.945 0.169 0.029 0.259 0.640 


60% 0.176 0.027 0.363 0.907 0.390 0.017 0.362 0.217 
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从 表 1， 我 们 可 以 得 到 如 下 结论 : 

(i) ” 随 着 缺失 概率 的 增 大 ， 利 用 ADEL 方法 和 利用 NAEL 方法 所 得 到 的 绝对 偏差 (Bias) 都 
逐渐 增 大 。 但 利用 ADEL 方法 所 得 偏差 的 增 大 速度 明显 小 于 利用 NAEL 方法 所 得 偏差 的 增 大 
速度 。 

(ii) ADEL 方 法 和 NAEL 方法 均 能 给 出 较 小 的 标准 误差 (SE)， 并 且 两 种 方法 所 得 的 标准 误 
差 对 缺失 概率 的 变化 都 不 是 太 敏感 。 另 外 ，NAEL 方法 给 出 的 标准 误差 略微 小 于 ADEL 方法 给 
出 的 标准 误差 ， 这 主要 是 由 于 ADEL 方法 中 的 调整 因子 含有 估计 量 缘故 。 

(ili) 尽管 利用 ADEL 方法 和 NAEL 方法 所 得 估计 的 置信 区 间 长 度 相差 不 大 ， 但 当 缺 失 概 
率 较 大 时 ，NAEL 估计 具有 较 大 的 绝对 偏差 ， 从 而 导致 其 覆盖 概率 是 非常 不 准确 的 。 

总 之 ， 从 模拟 结果 来 看 本 文 提出 的 ADEL 方法 具有 优良 的 有 限 样本 性 质 。 并 且 当 缺失 概率 
较 大 时 ，ADEL 方法 明显 优 于 NAEL 方法 。 


4 定理 的 证 明 


为 书写 方便 ， 下 文 用 c 表 示 正 常数 ， 每 次 出 现时 可 以 代表 不 同 的 值 。 在 证 明 本 文 的 主要 结 
果 之 前 ， 首 先 给 出 一 些 正则 化 条 件 。 

Cl: ”带宽 满足 h = Cn, HPO > 0 为 某 给 定常 数 。 

C2: RRK K) 是 对 称 的 概率 核 函 数 ， 且 有 


J AERA 


C3: 对 任 给 的 w € (0,1), f(u), (u), Y(u), o?(u) 以 及 0(w) 在 4 点 均 二 次 连续 可 微 ， 其 
中 f(-) 为 UU 的 密度 函数 。 
C4: 
sup E(e} |U: =u) < oo， sup E(X; |U: =u) < co， 
0<u<1 0<u<1 


ART uit, i=1,---,n, 7 = 1,… ,p， 其 中 Xi; EX, 的 第 r 个 分 量 。 
C5: MAH u, V(u) 是 正定 矩阵 。 
下 面 首先 给 出 几 个 引 理 。 
引 理 1 设 oaz :ani b1,b2,--- ,bn (by > b2 > … > bn) 为 两 实数 序列 。 令 


k 
Sk = 》 ai, 
i=1 


则 有 
| D aid; | <c max |b;| max |Sil. 
451 


1<i<n 1<i<n 
证 明 利用 Abel 不 等 式 ， 经 简单 计算 可 得 到 本 引 理 的 证 明 。 
引 理 2 We, i = 1,… ,n， 为 一 相互 独立 的 随机 变量 序列 ， 且 满足 BE(e;) = 0M Ele?) < 
C < OO, 则 有 





= Or(Vmnlogm). 


k 
max | > ei 
1<k<nl < i 

i= 
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进一步 ， S (ja, jz , jn) 为 (1,2,--- n) 的 任意 置换 ， 那么 有 





max 
1<k<n 


k 
Ye; | = Op(Vnlogn). 
i=l 


证 明 利用 Kolmogrov 不 等 式 ， 我 们 可 以 证 明 本 引 理 。 
并 且 
E\Yi|*° < 00, sup f lyl* f(a, y)dy < o, 


这 里 f(-,-) 是 (X,Y) 的 联合 密度 函数 。 设 天 (.) 为 一 有 界 支 撑 上 的 非 负 函 数 ， 并 满足 Lipschitz 条 
件 ， 那 么 有 


sup| PAL -BLOG — 2)Y;1}| = 0p( { EE), 


i 


FEF nh — oœ, 8< 1- s7 
证 明 证 明 见 Mack 和 Silvermanll6] 。 
引 理 4 设 条 件 Cl-C5 成 立 ， 则 有 


sup _||f(u) — ¥-*(u)®(u)|| = O,(C,), 
O<u<l 
sup jôu) = V1(u)B(wB 0(w))| = Op(Cn), 
O0<u<i 
Ht Alu) 和 9(w) 由 (5) 式 所 定义 ， 并 且 
cs 
证 明 结合 引 理 3， 利 用 类 似 文献 [9] 中 (7.1) 式 的 证 法 可 知 ， 对 we (0,1), 一致 有 


DIQ,Z = nA(u)f(u)®(u) ® (1,0)? {14+ 0,(C,)}, (10) 


T =n ; : 
Da Qu Du = A(u) f(u)U(u) = | 0 f s?K(s)ds 


{1+0,(C,)}, (11) 
其 中 A(w) = P(6; =1|/U, =u). BAB AW 的 定义 ， 经 简单 计算 可 得 
fi(u) = (Ip; 0p) (D3 Uu Du) D] OZ zz U~*(u)®(u) {1 十 O;(Cn)}, 


Xtu € (0,1) 一 致 成 立 。 利 用 类 似 的 方法 可 以 证 明 第 二 式 ， 这 就 完成 了 引 理 4 的 证 明 。 
引 理 5 设 条 件 C1-C5 成 立 ， 如 果 B 是 参数 真 值 ， 则 


Fa 8)  NO,B) 
i=l 


其 中 B= E{o?(U)[E(ZZ7) — 67(U)¥-1(U)B(U) }}. 
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证 明 结合 (7) 式 ， 经 简单 计算 可 得 
ĝi 


ô; 
市 (9) = 一 二 [2i - uU)" Xie + - 王 
A(U;) y Ai) 


2 [Z — u(U;)" Xi] XF [9(U;) — 9(Ui) + (Ui) 8] 


[ADD -ADD Xie 











[u(U:) -PC CAXE [0(Ui) — (Ui) + ACU) 6] 





+ 
AW) 


= Ja + Jiz + Jig + Jia. 
注意 到 在 条 件 C1, C2 下 有 


sup |A(u) — A(u)| = Op((nh)~2). 


再 利用 对 To E A(u) 点 进行 Taylor 展开 ， 并 经 简单 计算 可 得 








RL = ROE OD [Zi — (Ui) Xi) 
aan (A(Ui) — A(Ui)) [Zi — p(Ui)T Xi]es + op (1) 


le 1 oe 
= Ia + 一 I; 十 op(1). (12) 


TERE) EU) = 0, Var(Ila) = B 十 o(1)， 利 用 中 心 极限 定理 可 知 n-3 5 Ii N(0,B). F 
i=1 


iEn? X Tia 220: 令 
i=1 


k 0; 
b(U:) = A(Ui) - A(Ui), ai = NORLA — (Ui) Xi)e 

为 qx1 维 向 量 。 另外 ， 令 ais 为 ai 的 第 s 个 分 量 ， (b(U%;;), i= 1, pace n) 为 (b(U:), i= T; A , 1) 
的 某 一 置换 ， 使 得 b(Ui,) > b(U;,) 22 b(U;,,)> 对 应 的 (d,s, i=1,--- n) 记 为 (ati s, i= 
l,e ,nN)o 设 [2,s 为 Liz 的 第 s 个 分 量 ， 由 引 理 1 至 引 理 4 可 得 





1 
< — sup |b(u | a 
ee | Be > = 








ic lic 
[Ja ade = | DeWeese 


1 
= FO r(Cn)Op(Vmlogn) = op(1), 


即 n- Yo Lig 0。 再 由 (12) 式 知 n- $ Ja > N(0,B)。 
i=l tal 
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结合 {Xiei} = OM E{[Z; 一 py(0i)7Xi]X7} = 0 以 及 引 理 4， 用 类 似 证 明 (12) 式 的 方法 可 
以 证 明 n-3 S Jy 240, v = 2,3。 另 外 ， 由 引 理 4 得 
i=l 


[aE E a < Op vmOR) = op). 
i=l 


这 就 证 明了 本 引 理 。 
引 理 6 设 条 件 C1-C5 成 立 ， 如 果 是 参数 真 值 ， 则 


Lome 
= DAOA > B. 
i=1 
证 明 我 们 仍 使 用 引 理 5 的 记号 ， SJ? = Jiz + Jiz + Jia, BA 
1 DAON ele pA 十 L Sy PTH 1 DA + l SOT 
nag i S $ MEn iii We 7 es 


t=1 


Ay + Ag + Ag + Ag. 


Ill 


结合 (12) 式 知 


41 


二 2 Tal 十 z 2 Tata Pa 2 Tali +a Y lalh + op(1) 


i=1 
= Aj. + A12 + A13 + A14 + 0p(1). 


由 大 数 定律 可 得 An > B。 下 面 证 明 A - 0。 设 Aiz,rs 是 Ar 的 第 (7, s) 个 元 素 ，Iij; 是 Ly, 
了 = 1,2 的 第 r 个 分 量 。 利 用 Cauchy-Schwarz 不 等 式 可 得 


we (13) 
i=1 i=1 


由 引 理 5 知 
n`! Solar = Op(1)， nl > Ibs = op(1). 
i=1 i=1 


因此 ， 由 (13) 式 可 以 证 明 A > 0。 类 似 地 可 以 证 明 A), > 0, v= 3,4. Bt, Ay > E. 
利用 类 似 的 证 明 方法 ， 可 以 得 到 4。 = 0, v = 3,4， 这 就 证 完了 本 引 理 。 
定理 1 的 证 明 ”结合 引 理 5， 经 简单 计算 可 得 
max (BO) = p(n4). 
再 结合 引 理 5， 利 用 与 Qin F Lawless LA TEI AT Ml Al = Opn?) Ak, HRA 
FA Taylor 展开 可 得 


R(B) =2 > {rT H(B) — [AT 市 (6)]2/2 } + op(1). (14) 


i=1 
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再 对 (8) SLA Taylor 展开， 并 经 简单 计算 可 得 


2 [ATAB]? = > ni(B) + op(1), (15) 
A= { Vaal DAO) + op(n-4). (16) 
i=l i=1 


HH (14)-(16) 式 可 得 


其 中 F 
B =n Y AlO). 
i=1 


再 利用 引 理 5 和 引 理 6 就 证 完了 定理 1。 
定理 2 的 证 明 ”结合 (9) 式 ， 利 用 文献 [4 中 定理 2 类 似 的 证 法 可 得 


Ê- B= nD 市 (9) + 0,(n-2). (17) 


i=1 


利用 引 理 6 的 证 法 可 得 个 2 TT， 这 里 TT = E(ZZT) — B{BT(U)Y-1(U)@B(U)}。 结 合 引 
理 5 及 Slutsky 定理 ， 由 (17) 式 可 得 


Vn(6 — B) “+ N (0, TBT’). 
AARAA E = Br, ARETE. 
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Empirical Likelihood Inferences for Semiparametric 
Varying-coefficient Partially Linear Models with Missing 
Responses at Random 


ZHAO Pei-xin!, XUE Liu-gen? 


(1- Department of Mathematics, Hechi University, Yizhou, Guangxi 546300; 
2- College of Applied Sciences, Beijing University of Technology, Beijing 100124) 


Abstract: In this paper, we consider estimations of the semiparametric varying-coefficient partially 
linear models with missing responses at random. An adjusted empirical log-likelihood ratio function 
for parametric components is proposed. The Wilks’ phenomena is proved and the confidence regions 
for parametric components are constructed. The optimal convergence rate and the semiparametric 
efficiency bound of maximum empirical likelihood estimators are derived. Simulation results show that 
the adjusted empirical likelihood method outperforms the unadjusted empirical likelihood method. 
Keywords: varying-coefficient partially linear models; empirical likelihood; confidence regions; missing 
data 
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